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ABSTRACT 

Aims. We investigate the structure of the interstellar medium (ISM) and identify the location of possible embedded excitation sources 
from far-infrared (FIR) line and mid-infrared continuum emission maps. 

Methods. We carried out imaging spectroscopic observations of four giant Galactic star-forming regions with the Fourier Transform 
Spectrometer (FTS) onboard AKARI. We obtained [OIII] 88pm and [CII] 158/um line intensity maps of all the regions: G3.270- 
0.101, G333.6-0.2, NGC 3603, and M17. 

Results. For G3. 270-0. 101, we obtained high-spatial-resolution [OIII] 88/um line-emission maps and a FIR continuum map for the 
first time, which imply that [O III] 88 /urn emission identifies the excitation sources more clearly than the radio continuum emission. In 
G333. 6-0.2, we found a local [O III] 88 pm emission peak, which is indicative of an excitation source. This is supported by the 18yum 
continuum emission, which is considered to trace the hot dust distribution. For all regions, the [C II] 158 pm emission is distributed 
widely as suggested by previous observations of star-forming regions. 

Conclusions. We conclude that [O III] 88 //m emission traces the excitation sources more accurately than the radio continuum emis- 
sion, especially where there is a high density and/or column density gradient. The FIR spectroscopy provides a promising means of 
understanding the nature of star-forming regions. 

Key words. infrared:ISM — ISM: HII regions — ISM: lines and bands 

1 . Introduction l2001b iMizutani et ail 120021: lOkada et alJ I2003L 120061) . In the 

present paper, we present FIR imaging spectroscopic obser- 

_ . _ , . . vations with the Fourier Transform Spectrometer (FTS) on- 

Far-infrared (FIR) spectroscopy of star-forming regions pro- board AKART (Murakami et al 2007 Kawada et al 2008)i the 

vides us with a great deal of information about the inter- T ■ c , . • , . , , , ■ OAA( - 

„, , . ry , , , Japanese infrared astronomical satellite launched in 2006. The 

stellar medium (ISM) without being affected by the uncer- C rf, c • . • . f ^ T e A c 

_ , \ / . & „ T _ . , • ,, ,. FTS is the spectroscopic part of the Far-Infrared Surveyor 

tainty of the extinction correction. FIR forbidden lines are /C . TC , . , i h A n-7k c- ,u t ■ u-r* 

J . . . , . (FIS; Kawada et al. 2007). Since the spectroscopic capability 

less sensitive to the electron temperature than optica for- rf ^ M ultib and Imaging Ph otometer for Spitzer (MIPS) on- 

bidden lines in HII regions, and observations with balloon- , ■ c rrr, ^'TF^nA ■ * ■ * J* i ,l e 

, ., , 6 , . . board Spitzer (Werner et al. 2004) is restricted to wavelengths of 

borne and airborne telescopes have been used to determine _ nn t , , t - fT> u 1{ ,. ctc ct-c • 

. , . , F , , . , , < 97 urn at low spectral resolutions oiR — 15—25, the FIS-FTS is 

physical properties such as the electron density and elemen- • ,-. Tr) . ... . • e . , ■ Cln 

, , , , . ni , n „j rr- : — TirrxTTa k^a. a unique FIR spectroscopic instrument in space for studying FIR 

tal abundance jRubm etalJ [Mi | Simpson etal.| [1995] forbi g de n-Une emis sion after ISO/LWS ancJ before the Herschel 

The Long-Wavelength Spectrometer (LWS; Ckgg et a L || l99f|) § observatory (Pilbratt 2008). With two dimensional arrays 

onboard the Infrared Space Observatory (ISO; | Kessler et alJ t , FT c. PT c has thp Pnnnhi ,itv nf nht.inlna ^n^tr.l ™™ ™ m f 



unuuaiu u. c iniiaicu opa^c wu SC iva L uiy | xv Caa i C i ^ a H me FIS-FTS has the capability of obtaining spectral maps over a 

1226} provided continuous spectra of the wavelength range wMe ^ ^ a few ^ observations. 
43-197 pm, enabling us to investigate the physical conditions 
of the gas in the ionized and/or photodissociation regions 

(PDRs) by me asuring several emissio n lines and the contin- 2. Observations and data reduction 

uum emission^ (Martfn-HernandeziFiDSQOl see| Abergel et alJ 2J Observation parameters 
120051 and iPeeters et alJ 120051 for reviews). Since the LWS y 

was a single-beam instrument, raster-scan observations were The observations are summarized in Table Q] All observa- 

needed to investigate the spatial di stribution of emission spec- tions were carried out with the Astronomical Observation 

tra and physical properties (e.g. iRodriguez -Fernandez et alJ Template (AOT) of FIS03, which is a spectroscopic mode of 
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Table 1. The summary of the observations. 



Target 


Number of 


Observation ID 


Observation Date 


AOT 


Reset interval 


Resolution 




pointings 




(UT) 




[s] 




G3.270-0.101 


1 


5110044 


Mar. 20 2007 


FIS03 


0.25 


Full 


G333.6-0.2 


2 


5110052/53 


Mar. 3 & 4 2007 


FIS03 


0.1 


Full 


NGC3603 


3 


1403056/57/59 


Jul. 21,22,23 2007 


FIS03 


0.1 


Full 


M17 


1 


5110024 


Sep. 27 2006 


FIS03 


0.1 


Full 



the AKARI/FIS. The FIS-FTS adopts a Martin-Puplett-type in- 
terferometer (see iKawada et alj|2008l for details of the instru- 
ment). To investigate emission lines, we use the full-resolution 
mode with the spectral resolution of 0.19cm _I without apodiza- 
tion and 0.38 ctrT 1 with 'Hanning' apodization. We obtain 7 sets 
of forward and backward scans of the moving mirror with a to- 
tal ~ 11 minute integration time for each pointing observation. 
The FIS has a short wavelength array, SW, and long wavelength 
array, LW, each of which has narrow and wide bands. In the 
spectroscopic mode, only a wide band of each array is used. 
They are 3 x 20 and 3 x 15 pixel arrays that cover 85-130 cm -1 
and 60-88 cm" 1 , on pixel scales of 26.8"and 44.2", respectively 
(IKawada et al.ll2007l 120081, Murakami et al. in prep.). Some of 
the present observations were carried out as Director's Time ob- 
servations, and the others as the Mission Program of ISMGN 
(Interstellar dust and gas in various environments of out Galaxy 
and nearby galaxies: P.I. H. Kaneda). The field-of-view (FOV) of 
SW and LW is 1 .5' x 9.8' and 2.4' x 12.2', respectively. Because 
of the misalignment of the arrays, the FOVs of SW and LW are 
not completely matched with each other and cover slightly dif- 
ferent areas of the sky (IKawada et al.ll2007h . The point spread 
function (PSF) of the FIS-FTS can be described well by a com- 
bination of two Gaussian profiles, which are located along the 
major axis o f the detectors with t he separation of 14" in each de- 
tector array (IKawada et al.ll2008l) . The FWHM of each Gaussian 
is 39" and 44" along the minor and major axis, respectively, for 
SW, and 53" and 57" along the minor and major axis, respec- 
tively, for LW. The position accuracy of the FTS mode is less 
than a half of the pixel size (Murakami et al. in prep.). 

2.2. Targets 

To take advantage of the imaging spectroscopy of the FIS-FTS, 
we selected four active star-forming regions that are bright and 
extended enough to detect the interferogram with a high signal- 
to-noise ratio (S/N) over an area of ~ 10'. 

G3. 270-0. 101 is a giant HII region located 14.3 kpc 
from the Sun. It was observed in the radio and in the mid- 
infrared (MIR) with the Mid-course Space Experiment (MSX) 
dConti & Crowtherl |2004|) . The number of the Ly man contin- 
uum p hotons A^(Lyc) is estimated to be 3 x 1 50 s~ 1 . lOkada et al.l 
(2008) detected several MIR ionic fine-structure lines with the 
Infrared Spectrometer (IRS) onboard Spitzer in this region. We 
completed one pointing observation, covering the positions ob- 
served by the IRS. 

G333. 6-0.2 is a bright H II region at radio and infrared wave- 
lengths dGoss & ShayeJll970UBecklin et al.lll973UHvland et al.l 
1980; Fui ivoshi et al.ll 1 998L 1200 lL I2005L 120061) t hat is totally ob- 
scure d in the visible. It is located at 3.1 kpc dConti & Crowtherl 
2004) and likely to be ex cited by a cluster of O and B 
stars (iFuiivoshi et al.l [1998), behind which dust grains with 
a lar ge variety of tem p eratures emit strong infrared emis- 
sion dHvland et aT1ll980l) . IFuiivoshi etail (120061) discovered the 
complex structure of the central 10 arcsec region using high 



spatial resolut i on ra dio observations. Colg an et al.l (1 19931) and 
ISimpson et al.l (120041) detected several MIR and FIR forbidden 
lines and examined the electron density and the abu ndance ra- 
tio at the center and several northern positions. Okad a^etall 
(2008) performed MIR spectroscopy with the Spitzer/IRS as for 
G3. 270-0. 101, and detected 14 emission lines from metal ions 
and molecular hydrogen, indicating the coexistence of gas in var- 
ious excitation states. 

NGC 3603 is one of the most massive, optically visible H II 
regions in our Galaxy. It has an integrated Ha flux of L(Ha) 
~ 1.5 x 10 39 erg s _1 , equivalent to 20 05 V stars dKennicuttl 
11984 . It is located in the Carina spiral arm at a distance of 
about 7 kpc (Nurnberger et alj2002l and references within). The 
central cluster, HP 97950, contains 3 Wolf-Rayet (WR) stars 
(Mof fat et al.l Il994) and it is understood it experienced a sin- 
gle burst of star formation 1 Myr ago dStolte et alj|2004l) . It 
shows a remarkabl e similarity to th e dense core R136 in 30 
Dor in the LMC dMoffat et al.lll994l) . iNiirnberger eTail (120021) 
determined the properties and structure of the molecular gas 
using CS observations. The structure of the ISM and infrared 
sources were examined with M IR observations (Nurnberger 
120031: INiirnberger & Stankell2003l) . MIR spectroscopic observa- 
tions were performed using the Spitzer/IRS and spatial variations 
in both polycyclic aromatic h ydrocarbons (PAHs) and elem ental 
abundances were examined dLebouteiller et alj|2007l [2008). 

M17 is also one of the brightest giant HII regions in our 
Galaxy and has been extensively studied at various wavelengths. 



It con tains a dozen of O- or early B-type stars (Hanson et al 



1997) and the age of the central cluster is < 3 Myr dHanson et al 



1997UJiang efa l. 2002). iFelli et alj d 19841) presented radio maps 
of M17 and discussed the structures and nature of the ISM. 
The northern bar is clearly seen in optical emission, which 
is indicative of a small amount of foreground extinction, 
whereas the southern bar has no clear optical counterpart and 
is highly o b scured . Emission at [OIII] 88 pm was detected by 
War d et al.l d 19751) . and the electron density and nature of the 
ionizing sources were investigated based on observations o f 
the FIR [OIII] and [NIII] emission lines (lEmerv et al.lll983l) . 
The PC III 158 ^m emission is widely dist r ibuted in this region 
(IStutzki et al.lll988l:lMatsuhara et al.lll989h . lStutzki etafl fl988) 
investigated the dumpiness of the edge-on H II region/PDR in- 
terface using the spatial distribution of [C II] 158 pm and CO ro- 
tational line emissio n. The PDR structu res and morphology wer e 
also investigated bv lAndo etail. d2002l) and lPovich et al J (120071) . 



2.3. Data reduction 

We performed the data reduction of the time series data of the 
FIS-FTS using the official pipeline. Details of the processes and 
calibration are described in Murakami et al. (in prep.). The time 
series data are provided as integrated ramp signals. The interfer- 
ogram is derived from the differentiation of the time series data 
after correcting for the non-linearity of the ramp curve. We re- 
move glitches in the interferogram and apply the discrete Fourier 
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Fig. 1. (a) An example of the obtained full-resolution spectra in G333.6-0.2 region. Three emission lines are labeled. For the 
spectrum of 60-89 crrT 1 , obtained by LW detector, we apply the defringing process with fringe parameters of around [N II] 122 fim 
line and do not apply the apodization. For the spectrum of 88-135 cm 1 , obtained by SW detector, we apply the defringing process 
with fringe parameters around [O III] 88 yum line and apply the 'Hanning' apodization (see Appendix©, (b) The lower resolution 
version of (a) obtained by the Fourier transformation of the shorter range of the interferogram (see text). 



transform with the self-determined zero optical path difference 
to obtain the spectra. After the flux calibration, we have a spec- 
trum of each forward and backward scan of the driving mirror. 
We then take an average of 7 scans of the individual scanning 
directions. 

[OIII] 88/i/m is detected by the SW detector and [Nil] 
122yum and [Cll] 158 /mi are detected by the LW detector. The 
spectr a obtained by both detectors are affected by fringe pat- 
terns dKawada et al l 12008. Murakami et al. in prep.). The de- 
fringing process and the fitting function are optimized for each 
line. Details of the line-fitting process are described in Appendix 
lAl FigureQJ shows an example of a spectrum (an average of the 
forward and backward scans). Since the parameters of the de- 
fringing process are optimized for fringes around emission lines, 
the residual of the fringe pattern becomes large at wavenumbers 
far from the emission lines. After obtaining the line intensities 
for each pixel within the detectors, we combine the results from 
2 or 3 pointing observations of the same region (table [TJ and 
compile maps with a regular grid of half the size of the original 
pixel scale. Details of constructing line maps and estimating the 
uncertainties are described in Appendix [B] The final line maps 
are shown in Figs. [2] |7] [TT] and[14]for individual regions, which 
are discussed in the next section. 

The maps of the continuum emission at specific wavenum- 
bers are created by a similar procedure to that of the line maps. 
Instead of the defringed full-resolution spectra, we use low- 
resolution spectra, which are obtained by the Fourier transfor- 
mation of a narrower range of the optical path (+0.35 cm for 
SW and +0.4 cm for LW) to avoid the effect of fringes (Fig.QJi). 
The obtained continuum maps are also shown in Figs. [2l 171 ITTI 
and[H 
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Fig. 4. The calculated ratio of 120 crrr'/lOO crrT 1 versus the dust 
temperature with the dust emissivity index of -1, -1.5, and -2 
from the lower to upper curves. 



1 19991) is employed. The resolutions of these surveys are 45" and 
43", respectively, which are similar to those of the FIS-FTS. For 
the MIR continuum maps, we use 9 fim and 18 /mi dat a obtained 
by the All-Sky Survey of the Infrared C amera (IRC; Ona ka etaf] 
120071: llshihara et al]|2006l |2007l 120081) onboard AKARI. When 
either of the IRC bands data is saturated, we use the MSX A- 
band and E-band instead. The resolu tion of the IRC All-Sky 
Survey and MSX is 9.4" and 18.3" (lEgan et al.lll999l) . respec- 
tively. 



3. Results and discussion 

In the following, we discuss the structure and the nature of the 
ISM in each star-forming region from the line and continuum 
emission maps obtained by the present observations as well as 
radio and MIR continuum data. For the radio conti nuum map, we 
use N RAO VLA Sky Survey (NVSS) at 1.4 GHz (ICondon et al.l 
1998) if available, otherwise the 843 MHz contin uum by the 
Sydney University Molonglo Sky Survey (SUMSS; Boc k et al.1 



3.1. G3.270-0.101 

We present the highest spatial resolution maps [OIII] 88 /im 
and [Cll] 158 //m line emission (Figs. |2j;,d) and continuum 
(Figs.|2^,f) obtained for this region to date. Figures [S^-e show 
the correlation of both [OIII] 8 8 /mi and the FIR continuum 
emissions with the radio continuum from the NVSS at 1.4 GHz, 
and 9 fim and 1 8 fim emission detected by the IRC all-sky survey. 
The peak of the radio continuum is located at the peak of the FIR 



4 



Y. Okada et al.: Star-forming regions with FTS onboard AKARI 



£ 
o 

LU 
Q 

< -2 



(a) 



0.0 



4 2 -2 
A RA [min] 

[0111] 88um [10" 7 Wm 2 sr" 1 ] 

2.0 4.0 6.0 8.0 




2 0-2 
A RA [min] 

continuum at 120 cm" 1 [10 17 W m" 2 Hz" 1 sr" 1 ] 



1.0 



2.0 



3.0 



4.0 



5.0 



g 

I 

O 
LU 
Q 

< -2 



(e) 




£ 
o 

LU 
Q 

< -2 



(b) 




0.0 



4 2 -2 
A RA [min] 

[CM] 158|xm [10 7 Wm 2 sr" 1 ] 

2.0 4.0 



£. 
o 

LU 
Q 

< -2 



(d) 




6 4 2 -2 -4 
A RA [min] 

continuum at 70 cm" 1 [10 17 W m" 2 Hz" 1 sr" 1 ] 



o.o 



2.0 



4.0 



£. 
o 

LU 
Q 

< -2 



2 -2 -4 -6 6 4 2 -2 

A RA [min] A RA [min] 

continuum ratio 120 cm" 1 / 100 cm" 1 



(f) 




■E. 
O 
LU 
Q 

< -2 




2 0-2 
A RA [min] 



(9) 



-4 



-6 



Fig. 2. (a) The observed (gray) area of SW and (b) that of LW, (c) [OIII] 88 /an, and (d) [CII] 158 /mi intensity map and contin- 
uum maps at (e) 120 cm -1 and (f) 70 cm -1 and (g) the continuum ratio of 120cm /100cm in G3. 270-0. 101. The origin of the 
coordinates is RA = 17 h 53 ra 35 s .8 and Dec = -26°10'59" (J2000) dLockmarJ1989l) . 
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Fig. 3. The gray scale shows the [O III] 88 yi/m intensity ((a)-(c)), the continuum emission of 120 cm 1 (d), and the continuum ratio 
of nOcmVlOOcm 1 (e) forG3.270-0.101. In (a) the NVSS 1.4GHz intensity, in (b) the IRC 9 yum intensity, and in (c) IRC 18 yum 
intensity are overlaid in contours, respectively. In (d) and (e) the contours indicate the NVSS 1.4 GHz intensity. The contours of the 
NVSS are drawn from 0.02 Jy beam 1 in a 0.02 Jy beam 1 interval, and those of t he IRC 9 yum and 18 yum maps are on a linear scale 
of arbitrary units. Asterisks in (a) show the positions observed by the Spitzer/IRS (Okada et al. 2008). The origin of the coordinates 
is the same as Fig. [2] 



continuum emission (Fig.[3ji), but it has an offset from that of the 
[OIII] 88 /vm (Fig.[3k). This offset is larger than the uncertainty 
of the position, which is smaller than the grid size in the figures. 
Since it takes 35 . 1 2 eV to form 2+ , the [O III] 88 yum line emis- 
sion probably peaks close to the ionization source. However, the 
free-free continuum is only sensitive to the emission measure, 
not the ionization structure of the H II region. Therefore, it is 
likely that the ionization source is located at the [OIII] peak, 
and that the radio continuum (free-free) peak is just a local peak 



in the emission measure, not necessarily the location of an ion- 
ization source. The [OIII] 88 yum emission also shows a peak 
close to the radio peak, which indicates that some 2+ ions still 
exist, and the large emission measure compensates for the de- 
crease in the number of 2+ ions, although we cannot exclude 
another ionization source existing at the radio peak. This pic- 
ture is supported by several other observational results. First, the 
18 yum continuum emission, which is u nderstood to trace the hot 
dust distribution (llshihara et al.ll200"7l) . has a broad peak includ- 
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Fig. 5. The [S III] 33 fim intensity (plus points) and the intensity 
ratio of [Ne III] 36 //m/[S III] 33 fim (asterisk points) in G3.270- 
0.101 observed by the Spitzer/IRS. The latter is shifted slightly 
to the left for clarity. These four positions are shown in Fig. [3^. 
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Fig. 6. The spatial distribution of the [OIII] 88 pun emission in 
G3. 270-0.101. All the data of grids in Fig. |2}; are plotted, with 
the origin being the center of the brightest grids. 



Fig. 9. The distribution of the [OIII] 88 //m emission line along 
the solid line in Fig.|7};. 



optically thin conditions. Figure [5] shows that the peak of [S III] 
33 //m is located at the peak of the radio continuum, and the peak 
of the ratio [Ne III] 36 fim/[S III] 33 fim is similar to that of [O III] 
88 fim. 

In summary, the peak of the radio continuum is likely to be 
illuminated externally by the exciting source, and the position of 
the exciting source is probably at the peak of the [OIII] 88/L/m 
line intensity map. 

In Fig. [6] the spatial distribution of the [O III] 88 fim emission 
is shown, according to which the ex citation energy is investi - 
gated, assuming a distance is 14.3 kpc (ICon ti & Cro wtherl200' 
The Lyman photon flux of 3 x 10 50 s" 1 dConti & Crowtherll200 
is indicative of 23 07-type stars or their equivalent. Assuming 
an excitatio n source of 23 07-type stars, an electron density 
of 31 cm -1 (IConti & Crowfhei 20041) . an oxygen abundance of 
4.57 x 10~ 4 dAsplund et al.ll2005l) . and the ph otoionization cross - 
section and the recombination coefficient in lOsterbrockl (fl989). 
we derive the ++ Stromgren radius to be 4.6 pc, which is com- 
patible with the observed distribution of [OIII] 88 fim emission 
(Fig.[6]). As a check, we also compiled a model using a Salpeter 
initial mass function that matches the total Lyman continuum 
flux, and found that the ++ Stromgren sphere radius would then 
be about 4.7 pc, very close to our value derived above. 



ing the radio continuum peak, but also extends toward the south- 
west (Fig. [3};). Second, Fig. [3£ shows the continuum ratio of 
120cm _I /100cm _1 , which is a good tracer of a dust temperature 
in the range 20-50 K (Fig.[4]i. For example, for a ratio of 0.8, the 
dust temperature (T<j) is 30 K, and for a ratio of 1.0, T& ~ 40 K 
with the emissivity index of - 1 . The high temperature (high ra- 
tio) region is as widely distributed as the [OIII] and the radio 
continuum emission. Observations of MIR emission lines using 
the Spitzer/]RS have also shown a spatial distribution that sup- 
ports this hypothesis. Figure [5] shows the intensity of the [S III] 
33 fim line and the intensity ratio of [Ne III] 36 ^m/[S III] 33 fim 
at 4 observed positions indicated by asterisks in Fig. [3^. The 
line i ntensities are averages of the dithering pairs ( Okad a~et al.l 
2008). The ionization potential of S 2+ and Ne 2+ is 23.33 eV 
and 40.96 eV, respectively. Thus, the intensity ratio of [Ne III] 
36 fim/[S III] 33 //m indicates the ionization degree, if a constant 
elemental abundance within the region is assumed, whereas the 
[S III] 33 fim intensity depends on the electron density and/or the 
gas column density as the radio continuum emission does under 



3.2. G333.6-0.2 

The exciting cluster of G333.6-0.2 is located around the origin 
of Fig. [7] where both the SW and LW detectors are saturated. 
Figure[8]shows the correlation of both [O III] 88 fim and the con- 
tinuum emission ratio of 120cm _1 /100cm _1 with the 843 MHz 
radio continuum measured by SUMSS and the 9 fim and IE fim 
emission detected by the IRC all-sky survey. Globally, the line 
emission or continuum exhibit a large gradient in the observed 
area and become intense towards the exciting cluster. The [O III] 
88 fim emission has a local peak around (ARA, ADec) ~ (2'. 5, 
2') (Fig. [7};). Figure [9] shows the intensity distribution along the 
line shown in Fig. [7}:, which confirms that this local peak is 
real even taking into account the uncertainty in the line intensity 
(AppendixlBj. Is this local [O III] peak indicative of an excitation 
source? 

Although the radio continuum does not show a clear feature 
corresponding to the [O III] local peak (Fig. [8^), the 9 fim emis- 
sion appears as a wall-like structure that sandwiches the eastern 
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Fig. 7. (a) The observed (gray) and saturated (black) area of SW and (b) that of LW, (c) [OIII] 88 /an, (d) [CII] 158 /an, and (f) 
[Nil] 122 //m intensity map and continuum maps at (e) 120 cm -1 and (h) 70 cm -1 and (g) the continuum ratio of 120cm _I /100cm _1 
in G333.6-0.2. The o rigin of the coordinates is the radio main peak [RA = 16 h 22 m 09 s .5 and Dec = -50°05'58" (J2000)] 
(Fujiyoshi et al. Il2006h . The solid line in (c) is a cut, along which the distribution of [OIII] emission is shown in Fig. [9] 
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Fig. 8. The gray scale shows the [OIII] 88/mi intensity ((a)-(c)) and continuum ratio of ^Ocirr'/lOOcm" 1 (d) for G333.6-0.2. 
In (a), the NVSS 1.4GHz intensity, in (b), the IRC 9 /mi intensity, and in (c), the IRC 18 /mi intensity are overlaid in contours, 
respectively. In (d), the contours indicate the NVSS 1.4GHz intensity. The contours of SUMSS are drawn from 0.1 Jy beam" 1 to 1 
Jy beam" 1 in a 0.1 Jy beam" 1 interval and 2 Jy beam" 1 and 3 Jy beam" 1 . The contours of the IRC 9 /mi and 18/im maps are on a 
linear scale with arbitrary units. The origin of the coordinates is the same as Fig. [7] 
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Fig. 10. The observed [CII] 158/mi emission versus the [Nil] 
122 /mi emission in G333.6-0.2. The lines show the calculated 
intensity ratio assuming all emission arises from ionized gas, for 
two electron densities, 16 cm" 3 (upper) and 977 cm" 3 (lower), 
the electron tempera ture of 10 4 K, and the solar abundance by 
lAsplund et ai1d2005l) . 



and western side of the [O III] 88 //m peak with a hole at the 
local [OIII] peak (Fig. [8j3). Since the 9 /mi emission is domi- 
nated by PAH emission, this hole in the 9 fim emission is under- 
stood to represent a depletion of the PAHs near the excitation 
sources, and an enhancement at the PDR/H II region interface. 
The 18 /mi emission has a corresponding feature at the [O III] lo- 



cal peak (Fig. [8};), which is consistent with heating sources that 
produce hot dust. On the other hand, [CII] 158 /mi (Fig. |7}I), 
[Nil] 122 /mi (Fig.[7J), the radio continuum (Fig. [8}i), and the 
continuum emission ratio of 120cm" 1 /100cm" 1 (Figs. [7g, [8}l) 
do not exhibit a clear peak around the [O III] local peak. If the 
[OIII] peak corresponds to the excitation source as in G3.270- 
0.101, the dust temperature should be high and thus the ratio 
of 120cm" 1 /100cm" 1 emission should also peak there. On the 
other hand, [N II] and [C II] emission should not show a peak 
at the local [O III] peak since most N and C ions should be in 
higher ionization states. The radio continuum is sensitive to the 
density and/or column density (i.e, the emission measure) of the 
ionized gas, and therefore may not show a peak at the excitation 
source. However, the flux ratio of 120 cm"Vl00 cm" 1 is sensitive 
to the temperature in a narrow range, and thus the absence of the 
peak in the ratio at the [O III] peak may be understood by this re- 
gion having a much higher dust temperature than G3. 270-0. 101 . 
The 120cm" I /100cm" 1 range around the [OIII] peak is higher 
than that in G3.270-0.101 (Figs.|2g and|2k). Unfortunately, other 
continuum flux ratios available in our data will not improve our 
understanding because of the poor S/N in other spectral ranges 
and misalignment of SW and LW detector arrays. Hence, we 
surmise but do not prove that the [O III] emission indicates the 
position of the exciting source(s). The distribution of the 18 /mi 
emission supports this hypothesis. 

Assuming that there is an excitation source at the [O III] local 
peak, we estimate the ++ -Str6mgren radius in the same man- 
ner as for G3. 270-0 . 101. With an electron density of 200cm" 1 
(Sim pson et al.ll2004l) . it is estimated to be 0.9 pc, 0.45 pc, and 
0.2 pc for an 03, 05, and 07 star, respectively, which cor- 
respond to 62", 30", and 14", respectively for a distance of 
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Fig. 11. (a) The observed (gray) and saturated (black) area of SW and (b) that of LW, (c) [OIII] 88 /mi, and (d) [CII] 158/rni 
intensity map and continuum maps at (e) 120 cm and (f) 70 crrT 1 in NGC 3603. The origin of the coordinates is HD97950 [RA = 
ll h 15 m 07 s .3 and Dec = -61°15'38".52 (J2000)]. 



3.1 kpc. The width of the [OIII] 88 /mi emission emitting area 
in Fig. [9] is comparable with the FWHM of the SW in spectro- 
scopic mode, which is 39"-44", thus the observed width is an 
upper limit. Therefore, an 03 star is too early to be the excita- 
tion source at the local peak of [O III] 88 /mi. On the other hand, 
the total 88 /mi luminosity that each ++ -Str6mgren sphere pro- 
vides is 2 x 10~ 13 , 6 x 10~ 14 , and 1 x 10- 14 Wm 2 , respec- 
tively. The observed 88 /mi excess around the [OIII] local peak 
is ~ 8 x 10 _l4 Wm 2 . Therefore, with taking account of the un- 
certainty in the absolute flux calibration of the FIS-FTS (50%; 
Murakami et al. in prep.), a single 05 type star can consistently 



account for the radius of the Stromgren sphere and the total flux 
of the excess 88 //m emission. 

Figure [10] shows the correlation between the observed [C II] 
158 /mi and [Nil] 122 jum intensities. They are basically well 
correlated with each other. We calculate the intensity ratio in the 
ionized gas for an electron density of 1 6-977 cirT 3 , which is de- 
rived from [Nil] 1 22 #m/205 um at the central position and two 
northern positions (IColgan et al.lll993tlSirnpson et al.ll2004 . the 
electron t emperature of 1 4 K, and the solar abundance for C 
and N by Asp lund et al.l d2005l) . The results are overplotted in 
Fig. [10] Since the [CII] 158/mi/[NlI] 122 /mi ratio largely de- 
pends on the electron density, 16-91% of the observed [CII] 



10 



Y. Okada et al.: Star-forming regions with FTS onboard AKARI 



intensity comes from the ionized gas for the range of n e - 
16-977 cm 3 . The contribution of the ionized gas to the [CII] 
158 fim emission has been investigated in so me star-forming re- 
gions: 27% or 20+10 % in the Carina ne bula dOberst et all2 006; 
iMizutani et aTll2004l) . 20-70 % in S 171 dOkada et al.ll2003l) . and 
30-80% in the cr Sco region (lOkada et al.ll2006l) . Model calcu- 
lations indi^ate_fhat it is at least 10 % and up to 50-60% ac cord- 
ing to lAbell d2006l) . and ~ 10% for lKaufman et alj d2006l) . The 
present result is compatible with these previous studies, although 
it does not provide a strong constraint. 



3.3. NGC 3603 

Figures[TT][T2] and[T3lshow the distribution of [O III] 88 fim and 
[C II] 158 fim line emission and radio, MIR, and FIR continuum 
emission in NGC 3603. These figures indicate that the overall 
spatial distribution of these tracers is similar to each other; these 
line and continuum emissions have two peaks toward the south 
of the central cluster (the origin in figures) and at around (ARA, 
ADec) ~ (-3', 2'), although some pixels of the SW detectors 
are saturated aroun d the former posi t ions. The positions of the 
molecular clump in Nurnberger et al. (2002) indicated by the CS 
observations are shown in Figs.fTZh and[T3b. Their positions are 
close to the peaks of the radio, MIR, and FIR continuum emis- 
sion, indicating the coexistence of the molecular clumps and the 
ionized gas, and the warm dust. 

The central cluster is located at the origin of FigQj] Previous 
observations suggest that a low density bubble of about 50" 
or less in radius is formed ar ound it facing to two dense cores 
of sou th-east and south-west dBalick et al.lll980t IClavton|[T986l 
1990). Because of the spatial resolution of the present observa- 
tions, the bubble detection cannot be confirmed, but we speculate 
that Figs. [Tit and [T2l show that the radio, MIR, and FIR contin- 
uum emission correspond to the structures of these dense cores. 
Based on the same discussion we made for G3. 270-0. 101, [O III] 
88 fim should have a peak at around the excitation sources, but 
a peak is clearly not seen because of the saturation of the pixels 

(Kg. im). 

On the other hand, there is another peak at around (ARA, 
ADec) ~ (-3', 2') in the [OIII] map as well as the radio, MIR, 
and FIR continuum emiss i ons, w hich corresponds to MM6 and 
MM7 in Nurn berger et all d2002l) . It indicates that the excitation 
source is within the molecular cloud. Here only the [C II] 158 /im 
shows a slightly different distribution, which extends to the west- 
ern side (Fig. [TBI and is indicative of an extended PDR. 



3.4. M17 

Observations were made across the southern bar region and part 
of the northern bar region in M17. The derived line and con- 
tinuum maps are shown in Fig. [14] and the MIR continuum 
emission by MSX is shown in Fig. [15] The NVSS radio con- 
tinuum can be obtained only for part of the observed region of 
the FIS-FTS and not shown in the figures. Unfortunately, nearly 
half of the pixels of the SW detector are saturated in the south- 
ern bar region (Fig. [14k). Toward the east of the southern bar 
corresponds to an HII region and the [OIII] 88 fim emission is 
detected, whereas the area west of southern bar contains molec- 
ular clouds and lacks any [OIII] 88 fim detection. The contrast 
between the east and the west side of the southern bar is also 
clearly seen in the continuum ratio of ^Ocirr'/lOOcirr 1 (not 
shown in the figure), which is ~ 1 .2 in the east side and below 1 
in the west side, indicative of hotter dust on the east side (56 K 



and 40 K, respectively, with the emissivity index of -1 from 
Fig- |4j- At the eastern edge of the observed region of the [OIII] 
emission, an increase in the line intensity can be seen, which 
corresponds to the northern bar and is consistent with previous 
observations (lEmerv et al.lll983l) . The [CII] 158/mi line emis- 
sion is widespread in this region and peaks at the northern bar 
and the southern bar (Fig. [T4"]D. in a good agreement with pre - 
vious o bservations (Stutz ki et al.M 1988: Matsuha ra et al .1 1989). 
lAndo et all d2002l) pointed out that a double ridge of the north- 
ern bar is distinct in the MSX band- A and the 3.3 fim UIR band 
emission corresponds to the outer ridge, again consistent with 
band-A tracing PAHs, hence PDR interfaces. Figure [T5h shows 
that the [C II] 158 fim emission is also strong at the outer ridge, 
which indicates the presence of the PDR region. 

The FIR continuum emission shows a strong peak at the 
southern bar, but there is no significant peak at the northern bar 
(Figs. n"4"b.f). This indicates that the north ern bar has a much 
lower column density than the southern bar ( Elmegreen & Ladal 
[19761 IStutzki et al.l[l988t lAndo et al.ll200llWilson et al.ll2003l) 



4. Summary 

We have presented the results of imaging spectroscopic obser- 
vations of 4 giant Galactic star-forming regions with FIS-FTS 
onboard AKARI. We obtained [OIII] 88 fim and [CII] 158 fim 
line intensity maps of all the regions. In G3. 270-0. 101, we found 
that the [OIII] 88 fim emission peaks at the excitation source, 
whereas the radio continuum emission shows a peak at a region 
0.8'(3 pc) away. In G333.6-0.2, we found a local [OIII] 88 fim 
emission peak 3'(3 pc) away from the central cluster, which is 
likely to be an embedded 05 star. In NGC 3603, the distribution 
of the MIR, FIR, and radio continuum emission as well as the 
[O III] 88 fim emission are similar to each other. For all regions, 
[CII] 158 fim emission is widely distributed, as suggested by 
previous observations of star-forming regions. We have pointed 
out that the [OIII] 88 fim emission shows a stronger correla- 
tion with the 18/um continuum emission than 9 fim emission, 
which indicates that the 18/im emission originates in hot dust 
grains close to early-type stars, in contrast to the 9 fim emis- 
sion, which traces PAH emission, and PAHs are destroyed in 
the relatively hard UV environments in which 2+ forms. This 
work indicates that the [OIII] 88 fim emission traces the exci- 
tation sources more closely than the radio continuum emission, 
especially when there is a large density and/or column density 
gradient. The FIR spectroscopy provides a powerful tool for un- 
derstanding the nature of star-forming regions including the po- 
sition of the excitation sources. 
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Appendix A: Line-fitting procedure 

Full-resolution spectra of both the S W and LW detectors are af- 
fected by stron g channel fringes, w hich can be modeled by an 
Airy function ( Ka wada et alj 120081 Murakami et al. in prep.). 
The defringing process is included in the official pipeline, in 
which the spectrum is divided by an Airy function with the pa- 
rameter being optimized for each emission line (Murakami et al. 
in prep.). We apply them to the present study. Since the [OIII] 
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Fig. 12. The [O III] 88 yum intensity map (gray scale) with contours of (a) SUMSS at 843 MHz, (b) MSX band-A at 8 fim, and (c) 
band-E at 21 fim in NGC 3603. The contours of the SUMSS are drawn from 10~' Jy beam -1 to 10 0- Jy beam -1 in a 0.2 interval 
of index. The contours of MSX band-A is 5 x 10 -6 , 1 x 10 -5 , 2 x 10 -5 , 3 x 10 -5 , 4 x 10 -5 , 5 x 10 -5 , and 1 x 10 -4 Wm 2 sr"\ and 
those of MSX band-E is the sam e with the additional two contours of 2 x 10~ 4 and 3 x 10 4 Wm 2 sr~'. The asterisks in (a) show 
the molecular clumps reported by Niirnberger et al. (2002) as MM7, MM6, MM1, and MM2 from north to south. The origin of the 
coordinates is the same as in Fig.fTTI 



emission is strong and the line width is sufficiently smaller than 
the fringe period, we apply the 'Hanning' apodization to avoid 
strong tails around the emission line. After defringing, we derive 
the line intensities measured by Gaussian fits, assuming a linear 
baseline. For the [Nil] and [CII] emission, which are detected 
by the LW detector, we do not apply the apodization to improve 
the detectability since the line width is comparable to the fringe 
period. The fringe pattern around the [C II] emission is more ir- 
regular than that around [Nil], and we perform a defringing for 
the [Nil] emission but not for [CII] emission. The amplitude 
of the fringe pattern around [CII] line is less than ~ 30 %. We 
derive the [N II] and [C II] line intensities by fitting with a sine 
function, assuming a linear baseline after defringing and without 
defringing, respectively. We then take the average of the line in- 
tensities of the forward and backward scans. The uncertainty in 
the line intensities is estimated from the uncertainties in both the 
baseline and the differences between the intensity of the forward 
and backward scans. 

Appendix B: Constructing line maps 

To create a map of a regular grid from distorted detector array 
images, and to combine images of several pointings toward the 
same region, we use a Delaunay triangulation algorithm using 
an IDL routine, interpolating the values from nearby points. The 



coordinates of each pixel of the detector array are estimated from 
the attitude of the satellite with the alignment of the detector 
arrays, including the distortion in the optics. 

The uncertainty in the line intensity of each pointing obser- 
vation is estimated from the uncertainties in both the baseline 
and the differences between the intensity of the forward and 
backward scans as described in Appendix [Al We estimate the 
uncertainty in the constructed maps by simulations. We assume 
a Gaussian distribution with the width of the uncertainty in the 
line intensity as the probability function for each pixel value. For 
each trial, we randomly set the value at each pixel according to 
this probability function, and create a map of a regular grid. After 
a thousand trials, we have a certain distribution for the value at 
each grid of the constructed maps. We fit the distribution with 
a Gaussian profile for each grid and consider the width of the 
obtained Gaussian to present the uncertainty. 
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Fig. 14. (a) The observed (gray) and saturated (black) area of SW and (b) that of LW, (c) [O III] 88 fan, and (d) [C II] 158 fim intensity 
map and continuum maps at (e) 120cm" 1 and (f) 70cm" 1 inM17. The origin of the coordinates is NGC6618 [RA = 18 h 20 m 26 s and 
Dec = -16°10'36"(J2000)]. 
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Fig. 15. The [C II] 158 fim intensity map (gray scale) with contours of (a) MSX band-A at 8 fim and (b) band-E at 21 fim in M17. 
The contours of MSX band-A are drawn from 1(T 5 - 4 to 1(T 3 6 W m 2 sr 1 in a 0.2 interval in logscale, and those of MSX band-E is 
the same but to 1CT 3 W m 2 sr~' . The origin of the coordinates is the same as in Fig. [14] 



